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, studied and a geometrical interpretation given in terms of star-hypersurfaces in quan- 

tum phase space. Our approach deals with discrete and continuous spectra in a unified 
fashion and includes a systematic treatment of non-diagonal stargenfunctions. The 
formalism is used to obtain a complete formal solution of Wigner quantum mechanics 
in the Heisenberg picture and to write a general formula for the stargenfunctions of 
Hamiltonians quadratic in the phase space variables in arbitrary dimension. A variety 
of systems is then used to illustrate the former results. 
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1 Introduction 



The Wigner formulation of quantum mechanics [1] - [8] has become an important field of 
research. This is probably due to the fact that Wigner theory formulates quantum mechanics 
in terms of "classical-like" objects. Because of this, it is perceived by many as more intuitive 
than the standard operator formulation and has been used to successfully address a consid- 
erable number of problems in a variety of fields of research ranging from the semiclassical 
limit of quantum mechanics [5, 9, 10, 11, 12], quantum chaos [13, 14] and hybrid dynamics 
[15] to M-theory [16, 17, 18, 19]. 

Possibly, the main advantage of Wigner quantum mechanics is conceptual and stems 
from its remarkable relation with classical statistical mechanics. The theory is formulated 
in terms of the phase space Wigner function which for one-dimensional dynamical systems 
reads: 

F w (q, P ,t) = -^-W(\m >< m\) = ^-f dye-™ < q+^y\m >< W)\*-%V >> (1) 
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and is the Weyl transform of the density matrix \tp(t) >< ip(t)\. The Weyl map W : A(H) — > 
A(T*M) can also be used to obtain the phase space representation A(q,p) = W(A) of a 
generic operator A belonging to the algebra of linear operators A(J~i) acting on the Hilbert 
space Ti. Furthermore, using the Schrodinger equation and the Weyl map we can easily 
derive the dynamics of the Wigner function, [3]: 

F w (q,p;t) = [H(q,p),F w (q,p;t)] M , 

where H(q,p) = W(H) is the Weyl symbol of the quantum Hamiltonian H. The Moyal 
bracket, 

[A(q,p),B(q,p)] M = — (A(q,p) * B(q,p) - B(q,p) * A(q,p)) , 

and the star product * are both H- deformations of the algebraic structures of classical me- 
chanics (the standard product and the Poisson bracket, respectively) [8, 9]. They can be 
defined through the relations: 

W(A-B)=W(A)*W(B), wQ^[A,B^j = [W(A),W(B)] M , V7L, B e A, 

from which their explicit functional form follows immediately: 

the derivatives d and d acting on A and B, respectively. With these structures the Weyl map 
becomes an isomorphism between the Lie algebra of quantum operators (A(T~C),-, [ , ]) 
and the Lie algebra of phase space functionals (A(T*M), *, [ , }m)- 

This structure yields basic physical predictions through the average value, 

< A(q,p;t) >= J dqdpA{q,p)F w {q,p;t), 

and the marginal probability functionals, 

P(q(t)=q ) = fdpF w (q ,p;t) and P(p(t) = p ) = J dqF w (q,p ;t), (2) 

both displaying an impressive similarity with their analogues in classical mechanics. 

To produce more general predictions one has to introduce the stargenvalue equation, 
which is the Weyl transform of the eigenvalue equation [20, 21]. For a one-dimensional 
system this equation reads: 

W(A\a >< a\) = W(a\a >< a\) <^> A(q,p) * g a {q,p) = ag a (q,p), (3) 

where \a > is a general eigenstate of A with associated (non-degenerate) eigenvalue a and 
9a(q,p) = W(\a >< a\) is the stargenf unction associated to the same eigenvalue. The 
stargenfunctions g a (q,p) can then be used to obtain the probabilities that a measurement of 
a generic observable A(q,p) yield the value a: 

P(A(q,p;t)=a) = J dqdp g a (q,p)F w \q,p;t), 
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thus generalizing eq.(2). 

The solutions of the stargenvalue equation can literally be connected to all relevant 
structures of the Wigner theory and it is thus not surprising that this equation has been 
extensively studied in the past. Its properties were systematically described in [9, 16, 20, 
21, 22, 23] and solutions (both analytical and numerical) for several particular systems were 
presented in [16, 20, 23, 24]. There is also a vast literature on the subject of determining 
the semiclassical behavior of these stargenfunctions. In his seminal work [25] , Berry derived 
the general semiclassical form for the eigenstates of a generic observable, using the Wigner 
function associated with the corresponding WKB-wavefunctions. Ozorio de Almeida and 
Hannay [26] generalized this approach to higher-dimensional systems. More recently, the 
dynamics of these semiclassical states has been studied in depth, [27, 28]. In particular, the 
authors addressed the issue of whether it makes sense to replace the Moyal bracket by the 
Poisson bracket for these semiclassical states. More details about these issues can be found 
in [29]. 

In this paper we revisit the problem of obtaining exact solutions for a generic stargen- 
value equation. A general formal solution will be presented both for continuous and discrete 
spectra. Our approach includes a systematic treatment of non-diagonal stargenfunctions 
and higher-dimensional systems and leads to 1) a geometrical interpretation of the stargen- 
functions in terms of what will be named star hyper surfaces in quantum phase space, 2) 
a complete formal solution of Wigner quantum mechanics in the Heisenberg picture and 
3) a general expression for the stargenfunctions of Hamiltonians quadratic in the phase 
space variables in arbitrary dimension. In this context we also dwell on the possibility of 
obtaining integral representations for polynomials associated with certain solutions of the 
time-independent Schrodinger equation. 

This last result integrates and generalizes several results previously presented in the lit- 
erature, regarding solutions of the stargenvalue equation for specific quadratic Hamiltonians 
[9]. It also proves the efficiency of the formal solution as a calculation tool in a variety of 
systems. 

However, it is important to emphasize that the main advantage of our approach is con- 
ceptual. The formal solution and the Heisenberg picture provide universal rules and formulae 
that are applicable to very disparate situations. They played a key role in relating the De 
Broglie-Bohm and the Wigner formulations of quantum mechanics (cf.[30, 31]) and were 
systematically used to study the semiclassical limit of quantum mechanics in [25, 26]. We 
also believe they might be particularly relevant for the field of constrained dynamical sys- 
tems [32] (and in particular for general parameterized systems) where the characterization 
of the physical space of quantum states is the main conceptual and technical issue of any 
quantization program [33]. 

This paper is organized as follows: in section 2 we discuss the general nature of projectors 
in standard operator quantum mechanics and use the Weyl map to derive the formal solution 
of a general stargenvalue equation. We also prove that stargenfunctions are an ^-deformation 
of the Dirac distribution and introduce the concept of star hypersurface. In section 3 we 
formulate Wigner quantum mechanics in the Heisenberg picture and present its complete 
formal solution. In section 4 we derive the explicit form of the stargenfunctions of a generic 
quadratic Hamiltonian in arbitrary dimension. Sections 5,6 and 7 illustrate the previous 
results through three simple examples: one-dimensional Hamiltonians with discrete (section 
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5 - simple harmonic oscillator) and continuous (section 6 - linear potential) spectra and a two- 
dimensional example (section 7 - two-dimensional harmonic oscillator), where we compute 
the simultaneous stargenf unctions (both diagonal and non-diagonal) of the energy and the 
^-component of angular momentum. Finally, in section 8, we present our conclusions. 

2 Projectors and stargenfunctions 

The cases of continuous and discrete spectra will be considered separately. 



2.1 Continuous spectrum 

Let us start with a one-dimensional system and consider a hermitian operator A with non- 
degenerate continuous spectrum. Let \a > be the general eigenstate of A with associated 
eigenvalue a. The explicit form of the projector \a >< a\, which will be designated by 
A (A — a), is given by: 

\axa\ = A(A-a) = ( dke lk{A ~ a) . (4) 

Z7T J 

To prove this result let us introduce two general states: \<p > and \ip >. Using the represen- 
tation of A we have: 

< V|A(i - a)|0 > = da' da" < ^\a' >< a'\ J dke ik{A ~ a) \a" >< a"\<j> > 

= J da'da" < i>\a' > 5(a" - a') J dke tk{a "- a) < a"\<f) > 

= J da' < ip\a' > 5(a - a) < a'\<p >=< ip\a >< a\<p >, (5) 

from which the identity (4) follows immediately. 

Let us now consider a dynamical system of arbitrary (but finite) dimension. Let A be 
an observable with eigenstates \a,z > such that A\a,z >= a\a,z > and z is a array of 
continuous degeneracy indices. It is then also true that: 

j dz\a,z>< a ,z\ = ^J dke lk{A ~ a) = A(i-a). (6) 

In fact: 

— < VI / dke lk{A - a ^\<j) >= 
= J J da'dz' J J da"dz" < ij\a',z' >< a',z'\ J dke lk ^\a",z" >< a",z"\<p >= 
= j j da'dz' j J da"dz"5(a" - a) < rj}\a! , z >< a', z'\a", z" >< a", z"\(j> >= 
= Jdz' < ^|a, z'xa, z'\<p >, (7) 

which proves that the functional form of the projector A(A — a) is always given by eq.(4) 
independently of the dimension of the system. 
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We now consider a two dimensional system. Let B be a second observable also with 
continuous spectrum and such that [A, B] = 0. The set of simultaneous eigenvectors \a, b > 
(such that A\a, b >= a\a, b > and B\a, b >= b\a, b >) spans the Hilbert space of the system. 
The projector \a,b X a,b\ is given by: \a,b X a,b\ = A(A — a)A(B — b), a result that 
follows from: 

A(A - a)A(B -b) = Jdb' J da'\a, b' >< a, b'\a', b><a',b\ = 

= J db' J da'\a,b' ><a',b\5(a' - a)5(b - b') = \a,b >< a,b\. (8) 

The generalization to higher dimensions is straightforward: let {A i: i = l..n) be a com- 
plete set of commuting observables displaying continuous spectra. The set of eigenvectors 
\a±, .., cii, .., a n > (such that Aj|ai, .., a^, .., a n >= aj|ai, .., a^, .., a n >, Vi = l..n ) spans the 
Hilbert space of the system. The projector \a±, .., a^, .., a n X a±, .., ai, .., a n \ reads: 

\ai, .., a,i, .., a n >< ai, .., a h ..,a n \ = A(Ai - a 1 )..A(A i - Oj)...A(A n - a n ). (9) 

Notice that it is crucial to our approach that the set of operators {A^} forms a complete 
set of commuting observables. If this is not the case (like for the Hamiltonian of a classical 
chaotic system) then the former resolution of the problem is not valid. From now on we will 
always assume that {A} is a complete set of commuting observables. 

To proceed we address the problem of obtaining the explicit functional form of the non- 
diagonal elements \b >< a\. We start by considering the one-dimensional case and introduce 
the "translation" operator: T(X)\a >= \a + A >. If B is such that [A, B] = %h then 
T(A) = exp(—i/h\B), and thus: 

\b >< a\ = f{b - a )\axa\ = ±- [ dke~^-^^-) = J- / dk e -^ b - a ^ k(A ^ , 

2,71 J 2,71 J 

(10) 

where in the last step, we used the Baker-Campbell-Hausdorff formula. The operator 
\b X a\ will be denoted by A(A,b, a). It is trivial to check that A(A,b,a) satisfies 
A(A,b,a)A(A,c,d) = 5{a - c)A(A,b,d) and so / daA(A, b, a) A(A, c, d) = A(A,b,d). A 
small remark is in order: the operator B is not completely specified by the commutation 
relation [A, B\ = %h. For instance B and B + F(A) satisfy the same commutation relations 
for generic functional F. Different operators B lead to different specifications of the non- 
diagonal element \b X a\ (cf.(10)). If B is required to be hermitian all possible elements 
| b >< a\ (for different B and fixed b, a) are related by a phase factor which is an exclusive 
function of A = b — a. This freedom has no effect on the physical predictions of the theory. 
Moreover our future results are equally valid for all possible choices of B in eq.(10) (providing 
they satisfy the proper commutation relations) and thus from now on, we shall assume that 
B is fully specified. 

These results can be generalized to higher dimensions. We introduce the translation 
operators Tj such that: Tj(Aj)|ai, .., a^, .., a n >= \a±, .., aj + Aj, .., a n >, Vi = \..n. If {B iy % = 
l..n} is another complete set of mutually commuting observables satisfying = ih5ij 

then Tj(Aj) = exp(— i/hXiBi), and the general non-diagonal projector reads: 

|a 1 , .., a n X ai, .., a n \ = T±(a 1 — a\)...T n {a n — a n ) \a±, .., a n X a± : .., a n \ = 
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= exp {-i/K{(a{ - ai)A + ... + (a^ - a n )B n }} A(A l - ai )..A(A n - a n ) = 

= A(ii, a'i, ai)...A(i n , <4, a n ). (11) 

Moreover, it is easy to prove that the single projector, 

\a[ >< Oj| = y dai... J ddi-i j dai + \... J da n \a± : .., a' i: .., a„ >< ai, .., aj, .., a n | (12) 
is given by: 

|a- >< cij| = T(a^ — aj)|aj >< aj| = A(A,, a-, a^) (13) 

where |aj >< aj| is given by eq.(6). Therefore the non-diagonal single projector (12) is of 
the general form (10) in any dimension. 

Finally, let us consider the problem of determining the general solution of the stargen- 
value equation. We first consider the one-dimensional case. The most general stargenvalue 
equation for an arbitrary Weyl symbol A is written [20]: 

A*g ba = bg ba and g ba *A = ag ba (14) 

and is the Weyl transform of the corresponding eigenvalue equation in the density matrix 
formulation of quantum mechanics: A\ip >< <p\ = b\ip >< <f)\ and \ip >< <p\A = a\ip >< <f>\. 
The general solution of this last equation is given by the projector (10) and thus the Weyl 
transform of (10) is the general solution of eq.(14): 

A.(A(9,p),M) = - I dke-J (b - a)B(q > p) * ef = - f dk( r^-^^^-^) ^ 

2% J 2,% J 

(15) 

the star exponential e* being defined by: 

OO 1 

^ (w) = E-r^(9»pr» ( 16 ) 

where A(q,p)* n is the n-fold starproduct of A(q,p) and A = A(q,p) = W(A), B = B(q,p) = 
W(B). Furthermore, if a = b then eq.(15) yields the diagonal element which is of the form: 

A*(A(q,p),a,a) = A*(A(q,p) - a) = i- J dk^ A ^~ a \ (17) 

In section 2.3 we will prove that this object is formally a h- deformation of the Dirac delta 
function: A*(A(q,p) — a) = 5(A(q,p) — a) + 0(h), the full identity being valid for those 
observables satisfying A* n = A n . 

The generalization to n-dimensional systems is easily carried out if one uses the formal 
expression of the n-dimensional projector (given by eq.(ll)) as a starting point: 

Wflai, ..,a' n >< ai, ..,a n \) = A^(A 1 , a[, a t ) * .. * A*(A n , a' n , a n ) 

= A,(ii,..,A l ;flir- ) flI l ;ai,--,flr l ), (18) 

where A*(Aj, a[, ai) = W{A(Ai, a-, a^)} is the single stargenfunction also given by (15) this 
time with A = A(q±, .., q n ,pi, ■■,p n ) and B = B(qi, .., q n ,Pi, ■■,p n )- Moreover the notation 
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A*(Ai, .., A n ; a[, .., a' n ; a±, .., a„) was introduced to designate the most general n-dimensional 
stargenfunction. From the previous discussion one is led to the conclusion that: 

A*{A U ...,A 

A*(Ai, A n ;ai, ...,a^;ai, ...,a n ) * A = ajA*(Ai, A n ; a[, a' n ; ai, a n ), (19) 

an identity that is valid for all i = l..n and that can be checked explicitly by substitution of 
eqs.(15,18) in eq.(19). In particular, if = a[ then A*(Ai, .., A iy .., A n ; a[, .., ai, .., a' n ; a±, .., a^, .., a n ) 
is one of the a^-left and -right stargenfunctions of the observable A^ Furthermore, notice 
that the single stargenfunction A*(Ai, a\, a^ also satisfies the former stargenvalue equation 
(this time just for a single value of i). In fact, the relation between the single and the n- 
dimensional stargenfunctions is very appealing: on the one hand they are related by eq.(18) 
and on the other hand, from eq.(12) they also satisfy: 

A*(Ai, a[, ai) = J da 1 ... J da^i J da i+1 .. J da n A*(Ai, .., A h .., A n \ a ± , .., a-, .., a n \ a u .., a u .., a n ) 

(20) 

We conclude that in the context of Wigner quantum mechanics the n-dimensional stargen- 
functions can always be constructed from the single ones and therefore we shall henceforth 
focus on the one-dimensional case only. 



2.2 Discrete spectrum 

The case of discrete spectrum is slightly more involved. The first step will be to introduce a 
"continuous like notation" allowing for a formulation of the discrete spectrum case in terms 
of the continuous spectrum formalism. Using this notation the entire set of results of the 
last section can be easily translated to the discrete spectrum case. 

Let A be an observable with discrete spectrum and let {\a n >} form a complete orthonor- 
mal set of eigenstates of A with associated non-degenerate eigenvalues a n . We now introduce 
the "continuous like notation" by defining the continuous projector: 

\a><a\=J25(a-a n )\a n ><a n \, (21) 

n 

which is identically zero for all values of a that do not belong to the spectrum of A. The 
intention is to use the matrix elements \a >< a\ and the continuous spectrum formalism to 
reproduce the discrete spectrum results. We start by proving that {\a >< a\, a e 7Z} is a 
complete set of projectors. Let then \(p > and \tp > be two general states: 

< 0| da\a >< a\ip > = da^2 $( a ~ a n) < 4>\ a n >< OrelV' > 

J J n 

= < <t>\ a n >< On|^ > = < </# >, (22) 

n 

and thus / da\a >< a\ = 1. Moreover: 

| a' >< a'\a >< a\ = 5(a' — a n )5(a — a m )\a n >< a n \a m >< a m \ = 

n,m 



7 



= J25(a' - a n )5(a - a m )5 n>rn \a n >< a m \ 

n,m 

= 5(a — a') y^£(a — a n )\a n >< a n \ — S(a — a')\a >< a\, (23) 

n 

and thus, as expected \a >< a\ is a well denned projector. Finally, we consider the probability 
distribution resulting from using the continuous spectrum formalism. If \ip > is the state of 
the system then V(A = a) = tr(|^ >< tp\\a >< a\) = J2 n ^( a ~ a n)\ < i>\o-n > | 2 and thus: 



ra+e 

P(A = a) = limy da'V{A = a) 

ra+e 

= lim / da' 5{a — a n )\ < ip\a n > | 2 = 



(24) 



if a 7^ a n , Vn 



< ip\a n > | 2 if 3n : a = a n 



as it should. 

The primary result concerning the stargenfunctions of A is that the projector \a >< a\ 
(21) is also given by eq.(4), that is: 

\a >< a\ — 5 (a — a n )\a n >< a n \ = A(A — a). (25) 

n 

To see this explicitly we introduce two general states \tp > and \<p > and proceed as in (5): 
< <fi\A(A - a)\ip > = ^~Y1 < 4>\ a n>< a n \ I dkexp{ik(A- a)}\a m >< a m \ip >= 

^ n,m J 

= ^5(a m - a) < (j)\a n >< a n \a m >< a m \^ >= 

n,m 

= ^ZK a m - a) < <f>\a m >< a m \ii >=< (f>\a >< a\ip > . (26) 

m 

The straightforward corollary being that A (A — a) = if a ^ a n for all n. The generalization 
to higher dimensions follows exactly the same steps as in the continuous case and the discrete 
spectrum stargenfunctions also satisfy eqs. (6,8,9). 

The non-diagonal elements can also be easily obtained if one knows the explicit form of 
the translation operator T(A). Notice that in the discrete spectrum case this operator is not 
of the form used in eq.(10), given the fact that there is no operator B satisfying [A, B] = ih. 
For instance, for the harmonic oscillator and for A — H we have T(A = nwh) = a n and 
T(X ^ nwh) = 0, where a is the creation or the destruction operator and n G Z. 

In general, let a m < a n be two eigenvalues of A, and define A„ i?n = a n — a m . The 
translation operator T(A„ im ) is such that: 



A, T(\ nm ) ^n,rnT(\ n ,rn) ■ 



(27) 



We then have: 



A (f (A n , 



rn J 



n,m ) 



iflm An,m) T{\ n ^ rn ^\(l m > fl n T'(A njm ) |o m >, 
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which means that T(\ n ^ m )\a rn > is an eigenstate of A with eigenvalue a n . There is a wide 
range of operators that fall into this classification. In fact, let us consider the classical semi- 
simple Lie Algebras (A n , B n , C n , D n , Eq, ■ ■ •) and construct the Chevalley canonical form 
[34] of the generators: {E a , E- a , H a }. Here a G A represent the set of simple roots, and H a 
are the corresponding elements in the Cartan subalgebra. Then, for any given /-dimensional 
representation of the algebra, the operators E± a work as ladder operators in the sense of 
equation (27). The most celebrated example is the angular momentum associated with the 
su(2) Lie algebra, where E± a correspond to the ladder operators J± = J x ± iJ y , which 
implement translations in the spectrum of H a , or J z . 

Once we have found an operator T(X njTn ) satisfying eq.(27), we get: \a n >< a m \ oc 
T(\ n ,m)\ a m >< a m\- in the "continuous spectrum notation" the general non-diagonal ele- 
ment is then T(A„ irrt )A(A — a m ) and yields (using the Weyl map) the general stargenf unction 
of A(q,p), (let T(A n , m ) = W(f(X nim )): 

A.(A, a n , a rn ) = T(A n , m ) * i- f dkef A ^ a -\ (28) 

This is the most general formula for non-diagonal elements and it is applicable both to 
continuous as well as discrete spectra. 

A slightly more useful formula can nevertheless be derived, provided the steps \ n ,m are 
constant, i.e. |a n +i — a n \ = A, Vn. In that case, let T = T(A) be the translation operator. 
The spectrum is defined by a m = a n + (m — n)\, (to, n G /), where I is some set of integers. 
We can write: 

\a m >= \a n + (to - n)X >= P n ,mT m ~ n \a n >, (29) 
where (3 n>m is some normalization constant. Consequently, 

\a + n\ >< a\ = T n \a >< a\ = S(a — a,k)— \au + n\ >< at\- (30) 

kel yk,n+k 

The Wigner functions associated with the non-diagonal elements of the density matrix, can 
thus be read off from the previous formula. In this respect it will prove useful to compute the 
following product of exponentials: e aT ■ e tkA . Since the commutator T, A is proportional to 

T, we conclude that all commutators appearing in the Baker-Campbell-Hausdorff formula 
either vanish or are proportional to T. Let us then write: 

e af _ e ikA _ e ikA+a^(k)f ^g-Q 

where /i(k) is some function of k yet to be determined. If we expand the exponentials in the 
previous equation and equate powers of T and A, we conclude that: 

= (32) 



We then have: 



A(A - a) = — [ dk e af ■ e ik ^ = — [ dk e ^(A-a)+^t _ (33) 
2,71 J 2tt J 
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From eqs. (28,33), we obtain upon application of the Weyl-map: 

A,(A, a + nX, a) = (T*) n A.(A-a) = £ : \ef *A*(A- a) 



a=0 



(34) 



Notice that the previous formula is equally valid if the spectrum is continuous. We assume 
that the finite "jump" A = b — a takes place in N infinitesimal uniform steps of "length" e: 
a i+ i — di — e, i — 0, 1, • • • , N, with a = a, a^ = b and e = In that case, we have: 

f = exp(-^e.s) ~ (35) 
From the operator analog of eq.(34), we then have: 
A(A,b,a)=]im N ^ +00 ±fdk £re ih ( A - a ^ f \ a=o = lim N _ +00 ±£?e« f f dk e^"^ 

= lim^+oo^ (l - i iJ ^B) N Jdk e^-") = ± e ~^ b ' a ^ J dk e lk ^ = 

= ±jdk e -kiP-)B + ik{A-^) 

(36) 

where we used the fact that: 

/ x\ N 

lim 1 + - =e x . (37) 

2.3 Basic properties of the stargenfunctions 

In this section we summarize several properties of the stargenfunctions: 

1) Let A(q,p) be a real symbol with non-degenerate spectrum. Then both the observables 
and the Wigner function can be expanded in terms of the functionals A*(A, b, a). For the 
Wigner function the expansion reads: 

F w (q,p) = J dadb (/ dq'dp'F W (q',p')AM(q',P'),a,b^ A.(A(q,p),b, a), (38) 

and equally for a general observable: 

X(q,p) = J dadb (| dq'dp'X(q' ' ,p')A*(A(q' ,p'), a,bf) A*(A(q,p),b,a). (39) 

Furthermore, if X(q,p) = A(q,p) then eq.(39) reduces to: 

A(q,p) = J da aA*(A(q,p) - a), 

this being the inverse formula of eq.(17). 

2) The probability functional for a general observable is given by: 

P(A(q,p) = a)=J dqdpF w (q,p)A*(A( q ,p) - a), 
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and fully copies the analogous object of classical statistical mechanics. Moreover, the distri- 
bution A*(A(q,p) — a) satisfies: 

/ dqdpA(q, p)A*(A(q, p) - a) 

$ dqdpA*{A{q,p) - a) °" ^ } 

for all a belonging to the spectrum of A(q,p). These two properties suggest that the prob- 
ability of finding the observable A with the value a is given (just like in classical statistical 
mechanics) by the integration of the Wigner distribution function F w over the phase space 
hypersurface A(q,p) = a. In fact one has to be more careful: the distribution A* (A — a) 
does not in general identify the hypersurface A = a, due to the non-local nature of the star 
product: A* (A — a) is a "star delta function" which in general also assumes non zero values 
in phase space points not belonging to the hypersurface A = a. What can be said is that 
the star delta function A*(A — a) identifies the hypersurface A = a in the star phase space 
or, in other words, that it identifies the star hypersurface A = a. 

3) We now prove that A*(A) is a ^-deformation of the Dirac delta function. From eq.(16) 
and the definition of the star product it follows that: 

ef =V^f- V — AJ m >A...AJ m ^A, (41) 

d d d d 



where J = [Jf^-f^ — J^J- Let now mi + ... + m n _i = s. From eqs. (17,41) it is clear that 
A* (A) can be cast as a power series in h: 

MA)=±(^)'jdk±^l £ -J-,^..^ (42) 

s=0 V Z / J n=0 IL - m 1 +..+m n - 1 =s ''*l--''*n-l ■ 

and is also trivial to check that the zero order term (s = 0) of the previous expression is just 
S(A). Hence, we conclude that A* (A) is an ^-deformation of 5(A). 

4) Let us calculate the explicit form of the former expansion up to the third order in h. 
From the definition of A* n (eq.(16)) we get: 

A *n = yltll) y — -AJ mi A...AJ m - 1 A 

f=o V 2 / mi+..^.,=. mi!..m n _i! 

= An + l(j) e ^ n ~ V An ~ 2 + I (y ) 2 ^-l)(^-2)A«- 3 + O(^)(43) 

where we used the notation: 0i = \AJ 2 A = f^0- - (|^|) 2 and ©2 = \{AJ 2 A 2 - 
2A(AJ 2 A)} = f# (^Y - 2|^-f^ + f# f^) 2 - It follows that: 

v / j H ^ \ r)r> aaov aa av av z \ oa 



dq 2 \ dp ) dqdp dq dp dp 2 \ dq J 

e - = 1+ i(|) 2 (tt)2ei + I(f) 2 (^e 2 



e lkA + 0(h 4 ), (44) 



and therefore: 



MA) = S(A) - ^"(A) - ^e 2 S"'(A) + 0(h% (45) 
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Both eq.(42) and eq.(45) are in perfect agreement with previous derivations of the semiclas- 
sical form of a generic stargenfunction [7]. It is however important to remark that these 
semiclassical expansions should be taken with caution. The expansion is useful if it trun- 
cates at some order. If, however, this is not the case then the expansion may not converge 
order-by-order. There are few known cases which can be dealt with appropriately (see for 
instance [25, 26, 29]). 

3 Wigner quantum mechanics in the Heisenberg pic- 
ture. 

The results of the previous sections lead to a complete formal solution of Wigner quantum 
mechanics in the Heisenberg picture. In this scheme the time evolution of a general observable 
A(q,p) is given by the equation of motion: 

^- t A(q,p;t) = [A(q,p;t),H(q,p)] M , (46) 
which displays the formal solution: 

+oo « 

A(q,p;t) = ]T -[...[A(q,p;0),H(q,p)} M ...,H(q,p)] M = U^)' 1 * A(q,p; 0) * U(t), (47) 

n=0 n - 

where U(t) = e ~ itH( - q ' p ^ h [ s the time propagator. Moreover, the general stargenfunction of 
A(q,P',t) is given by: 

A,(A(q,p;t)-a) = ^- J dkef A b,P*)-«) = j/^-i * A,(A(q,p; 0) - a) * U(t), (48) 
and thus it equally satisfies the time evolution equation: 

^A*(A(q,p;t)-a) = [A*(A(q,p; t) - a), H(q,p)] M . (49) 

The previous equation together with the relation 

A(q,p;t) = J da aA*(A(q,p;t) - a) 

lead to the conclusion that the stargenf unctions A*(A(q,p;t) — a) encapsulate the entire 
information concerning the time evolution of the system. In particular, the probability that 
a measurement of A(q,p; t) at time t yield the value a is given by: 

P(A(q,p;t) = a) = J dqdpF w (q,p)A if (A(q,p;t) - a), (50) 

and satisfies the following suggestive formula: 

P(A(q,p;t) = a) =< A*{A(q,p; t) - a) > . 
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Consequently: 

^P(A(q,p;t) = a) = J dqdpF w \q,p)[A t (A(q,p;t) - a), H(q,p)] M = 

= < [A*(A(q,p;t)-a),H(q,p)] M > . (51) 

Finally, notice that the probability predictions can be easily connected with the correspond- 
ing formula in the Schrodinger picture if one notices that: 

J dqdpF w (q,p) (Uit)- 1 * A*(A(q,p; 0) - a) * U(t)) 
= Jdqdp(u(t)*F w (q,p)*U(t)- 1 )AM(q,P;0)-a). (52) 

4 Solution for quadratic phase space functionals 

In this section we will derive a general expression for the stargenf unctions associated with the 
Weyl symbol of a generic quadratic operator. We will address the particular case when the 
Hessian matrix of such an operator is proportional to a symplectic matrix. In this instance 
we will analyze thoroughly continuous and discrete spectra. All these formulae are applicable 
both to diagonal and non-diagonal stargenf unctions and they generalize the results of ref. [9]. 

The calculation of the *-exponential is the most crucial step. Following a technique de- 
veloped in ref. [9], we can derive the *-exponential for any polynomial of degree 2 in arbitrary 
dimension. To be more precise, let: 

A = z T Az + b T z, (53) 

where A is a symmetric, non-singular, 2N x 2N matrix, the superscript T stands for matrix 
transposition, and: 

z T = (Pir ■ • ,PN,qir ■ ■ ,Qn) , b T = (&!,•• -,b N , b N+1 ,- ■ ■ ,b 2N ) ■ (54) 

The restriction that the matrix A be non-singular is not indispensable as we shall see in the 
example of the linear potential. We then prove the following theorem. 

Theorem: The non-commutative exponential is given by: 

el A = (detcosfir 1/2 exp j-i (V + \b T A' 1 ^ JS2 1 JtanBJS A {z + ^A _1 bj - jb T A^ , 

(55) 

where 

B = h/3JS A JS^J, (56) 

J is the 2N x 2N symplectic form: 

'-(:-;) <") 
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and Sa is a 2N x 2iV matrix which satisfies: 

S^S A = A. (58) 
The corresponding *-genfunction then reads: 

A*(A-a) = ^Jdk [det cosB(k)}- 1/2 x 

x exp {-I (z T + I^A- 1 ) JS^Jtan^A;) JSU (* + l^ft) -ik(a+ \b T A~ l b)) , 
where -B(fc) is given by (56) with (3 replaced by ik. 



Proof: Let us first calculate e^ v with V = z T Az. According to (16) we have: 



+oo on 



?1=0 



where: 



Q n+ i = V*Q n 



V + —{JAJ)ijdidj 



(59) 



(60) 



In this expression, sum over repeated indices is understood and {<%} = (jjj^, • • • , • • • , ^ 

It is easy to check that the imaginary terms in the *-product (i.e. terms proportional to odd 
powers of H) do not contribute to eq.(60). 

Let us now consider a function ((3, z) which satisfies: 



|| = V<P + f {JAJ) i3 d t d j( p, 0(0, z) = 1. 



(61) 



If we expand ((3, z) in powers of (3, ((3, z) = J2n=o ^T-^n( z ), and substitute in eq.(61), we 
then get the following recursive relations: 



K n+ i(z) 



ti 2 

V+—(JAJ) ij d i d j 



K n (z), K (z) = 1. 



(62) 



Upon comparison with (60), we conclude that tt n (z) = K n (z), and hence: (/3, z) = e^ v . 
To obtain the non- commutative exponential, we thus have to solve (61). We consider the 
following Ansatz: 



((3, z) — M (13) exp z 1 A ((3) z 



(63) 



where M (f3) is a scalar function and A (f3) is a symmetric 2N x 2N matrix-valued function. 
If we substitute (63) in (61), we get: 



% = %Tr (JAJA) , 
A' = A + h 2 AJAJA, 



(64) 



where M' (13) = <gf and A' ((3) = Notice that: 



d/3 m 



J- 1 = J T = -J. 



(65) 
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If we define the matrix B as in (56), it then follows that B and tanS are antisymmetric 
matrices. Likewise we define the matrix A as: 



A = --JS^JtanBJSA- 

h A 



A T = --S^JttmBJ (S^Y J = ^rJS^J J S A J S T A J tan B J (s^ J S A X J 



This matrix is symmetric. Indeed, we have from (65): 

h~ A ~ " ^ A ) " h 

However, since B = h(3JSAJS A J, we conclude that: 

J S A JS T A J tan B J [S A l ) T JS A l J = tan BJS A JS A J 2 (S^Y JS A l J = - tan 5. 
Eq.(67) then reads: 

A T = —JS^JtemBJSA = A. 



(66) 

JS A . (67) 

(68) 
(69) 



From (66) we have: 

A' = - JS^-PSaJS^J (cos 2 B) ~* JS A = S T A J (l + tan 2 B) JS a = A- S^J tan 2 BJS A - 

(70) 

On the other hand we have: 

h 2 A T JAJA = S^JtanBJ^S^Y J 2 SlS a J 2 S^ 1 J tan BJS a = - S T A J tan 2 B J S A , (71) 
and we conclude that (66) is indeed a solution of (64). It remains to prove that: 

M{(3) = (det cos B)~ 1/2 . 

For a generic matrix G, we have: 



d_ 

dp 



det G (P) = det G (f3) Tr 



(72) 
(73) 



It then follows that: 

{gM(0) = ^ (det cos 5 (/?))^ 2 



-±M 3 (det cos B) Tr [(cos B)' 1 ± cos B ((3) 
= \MTr [tan B JS A JSl J 



Moreover we have: 



(74) 



2" v ' 2 V ' 2 V ' ' J 

where we used the cyclicity of the trace. Comparing (74) and (75), we conclude that (72) is 
indeed a solution of (64). 

Finally, we prove eq.(55). Let us define v — z + ^A _1 b. We then have: T>{v) = v T Av = 
A(z) + \b 7 A~ l B. Notice that the translation z — > v does not affect the *-product. Conse- 
quently: 



h 2 h h 

—Tr(JAJA) = --Tr (.J AJ 2 S A l J tan B J S A ) = -Tr (tan BJSaJS^j) , 



J)A(z) 
-*(z) 



'*(*) 



v T Av _ ^ h T A -l h 



, (76) 



where we have used the notation *(z), *(v) to specify the variables with respect to which we 
compute the *-products. From the expressions for A and v, we recover eq.(55). n 
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4.1 Determining the matrix Sa 



A few comments are in order concerning the existence and uniqueness of the matrix Sa 
defined by eq.(58). First of all, notice that Sa is not unique. In fact it is defined modulo some 
orthogonal matrix: if O is a 27V x 2N matrix such that 00 T = T = I, then if A = S^S A , 
we also have A = S'aS' a where S' A = OSa- On the other hand, any quadratic form Q(z) 
can be algebraically diagonalized, which means that if Q(z) = J2i<i<j<2N a ij z i z ji then there 
exists a set of variables Vi such that Q(z) = Q'(v) = Y%=\ ^i v h where {Aj, % — 1, • • • , 2 A?"} 
is a set of complex numbers (possibly, not all distinct). This can be expressed in matrix 
language as follows: 

v = P~ l z, Q(z) = z T Az = v T Dv = Q'(v), 

where: 

D = diag (X u ■ ■ ■ , \ 2N ) = P T AP. 

Notice that in the previous expressions, in general P T ^ P^ 1 and P is a complex matrix. 
The exception is when A is real, in which case P is a real orthogonal matrix. 

To proceed we define the matrix = diag (\/Ai, • • • , v^2~/v) • The matrix Sa then reads: 

S A = D^P- 1 . 

In some cases, when diagonalizing the matrix in this fashion is too complicated, we can alter- 
natively determine Sa by a straightforward generalization of Cholesky's method to complex 
matrices. In that case, we look for a decomposition of the form Sa = U, where U is an upper 
triangular complex matrix. Its elements can be computed iteratively according to: 



Ui 



— (an ^k=i u ki) i 



2N, j = i + l,---,2N. 



— — Y^'" 1 



To specify when this method is applicable to a complex matrix A let us define the set of 2N 
matrices Ai (i — 1, • • • , 2N): 



a ll a 12 

d\2 0,22 



( 



A,= 



an O12 Ol3 



Ai = (au), ,A 2 = 
The method is then valid whenever det(Aj) ^ 0, Vi = 1, • • • , 2N. 



a\ 2 022 «23 
\ Ol3 a 23 a 3 3 / 



4.2 A particular case 

The formula (55) can be drastically simplified in the following particular case. Suppose the 
matrix Sa is proportional to a symplectic matrix, i.e.: 

S A JS T A = aJ, (77) 
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where a is some nonzero constant. If a matrix is symplectic, then so is its transpose. It then 
follows that: 

S T A JS A = aJ. (78) 

Since S A S A = A, we get: 

AJA = S^SaJS^Sa = aS T A JS A = a 2 J. (79) 

Moreover, 

S A JS t a = aJ S A JS T A J = -al, (80) 
where / is the 2N x 2N identity matrix. Consequently, 

Sa 1 = --JS T A J, A' 1 = -\jAJ. (81) 
a a 1 

From eq.(56), we have: 

B = h{3JS A JS T A J = -ha[3 J. (82) 

And so: 

/ B 2 " = (iha(3) 2n I 

\ B 2 ^ 1 = i (tha(3) 2n+1 J, 

from which follows: 

sin B = i sin (ihaf3) J, cos B = cos (ihctfi) I, tan B — % tan (iha/3) J. (84) 
We then have: 

A = -^JSa 1 Jt&nBJS A = tan (iha(3) J 2 S T A J A S A = — ^ tan (iha(3) A, (85) 
n an an 

and 

M = (det cos By 1/2 = [det cos {ifia(3) I}~ 1/2 = [cos (ika(3)}~ N . (86) 
Equation (55) thus reads: 

ef 4 = [cos (ihaP)] ~ N exp { -4r tan (iha[3) + -J—6 T A _1 6 [tan - ) , (87) 

where A~ l is given by (81). 

How interesting is this specific case? To begin with, it applies to the harmonic oscillator 
in arbitrary dimension. Furthermore it is valid for any one-dimensional system. Indeed, 
a 2 x 2 symplectic matrix is any matrix of determinant 1. So any nonsingular matrix is 
proportional to a symplectic matrix. On the other hand, even if the 2x2 matrix is singular 
(which entails a = 0), we can still apply the previous results if we replace the vanishing 
eigenvalue by some arbitrary infinitesimal parameter. This will be our approach in the case 
of the linear potential (cf. section 6). 

Let us then consider some examples. The ensuing analysis will slightly generalize the 
results of ref. [9]. We assume that eqs. (77,78) hold for a ^ and that 6 = 0. We now prove 
the following lemma: 
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Lemma: If A = S a Sa is hermitean, then a is either real or purely imaginary. 
Proof: If A is hermitean then it follows from (81) 1 : 



SaSa — S A S A 



qT ct C* Q— 1 



a 



- S A S A J S A'J' 



From this we infer that: 

Sa = —-JSaJS a S% 

< S* A = — J S* A J S A S A, 
k S\ = —-^:S A S A JS A J. 

If we hermitean conjugate eq.(77) and use the previous equation , we get: 

S* A JS A = «7^(^) 2 JS* A JSlS A JS T A S A JS\j = a* J 
j^JS* A JS^JS A JS A J = a* J (^) 2 JS* A J 3 S A J = a* J 



(88) 



(89) 



(90) 



^-J = a*J 



a" 



(«*)' 



±a*. D 



If a = 0, then we have to regularize the matrix A as we shall see in the example of the 
linear potential. 

If a is real (and nonzero), then the *-exponential in eq.(87) yields for (3 = ik [9]: 



ikA 



[cos (hah)] N exp 



^ tan (ftccA;) 



(91) 



For fixed A ^ 0, this can be viewed as an analytical distribution in the variable k e 
C\|(2n+ l)^, ^ £ Being a periodic distribution (with period ||), it can be ex- 
panded in a Fourier series: 

+00 

— sfkA \ ' /'• iix\hl; 



where the Fourier modes are given by: 

ah ("ik 



inahk 



(92) 



(93) 



If we make the substitution z = e tahk , we get: 



2^ 

C„ — - — : 



2m Jc 



dz z- n+N (z 2 + exp 



N 



A z 2 -l 



ah \ z 2 + 1 



(94) 



Notice that the integral in eq.(93) is improper because of the singularities at k — i^n - This 
can be dealt with, if we choose C to be some contour in the complex z-plane containing 
the origin in its interior, but not the singularities z = ±i. Cauchy's theorem shows that C n 
vanishes for n < N and for n even. A straightforward calculation yields: 



a 



2n+N+l 



^(-DV-Aif- (§*) , 



n 



0,1,2, 



(95) 



1 The superscripts * and f denote complex conjugation and hermitian conjugation, respectively. 
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where L%(x) are the Laguerre polynomials. This can be proven by using the following 
representation of L"(x)[35]: 



k=0 



An interesting fact of this approach is that, as a by-product, we obtain an integral represen- 
tation for the Laguerre polynomials. Indeed if we equate eqs.(94) and (95), we get: 

K(x) = jf dz z- 2n ~\z 2 + l)— 1 exp , n, a = 0, 1, 2, • • • (97) 

A remark concerning the normalization is now in order. In the previous equations the 
coefficients C 2 n+N+i are not the normalized Wigner functions. In fact we have: 

J d N q J d N p C 2n+N+l (q, p) = (2nhf , n = 0, 1, 2, • • • (98) 

Note that this is consistent with the fact that C 2n +N+i = W (\n >< n\) and that the Wigner 
function is F^ n = ^ 2 J^ N W (\n >< n\). Consequently: 

« = ^-ALf-(^),„ = 0,l,2... (99) 
From the previous analysis, eq.(17) reads: 

A, (A - a) = — / dk <j)(k)e- tka = C2n +N+ iS (a - a n ) , a n = ha(2n + l). (100) 
27r J „=o 

If, on the other hand, a = 27 (7 e 9? \ {0}), we have: 



A* (A — a) = ^ J dk [cosh (fryfc)] N exp ^ tanh (^,7/c) — -i/ca 

_ 2 iv-2 r(f + ^)r(f-^) A F (N, Ja^. N . 2iA) 



(101) 



where 1-F1 is the confluent hypergeometric function [35]. 

These results are a slight generalization of those in ref[9], where Bayen et al. derived the 
expressions for continuous and discrete spectra, when the matrix A takes the form: 

A _ I 0/2x2 ^2x2 \ 

\ bl 2x2 c/ 2 x2 / ' 



for arbitrary real constants a, b, a This corresponds to a — V 'ac — b 2 in eq.(77), which may 
be real or purely imaginary. Obviously eq.(77) admits more solutions for A than the previous 
ones. 
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5 The one-dimensional harmonic oscillator 

The aim of this section is to determine the diagonal and non-diagonal Wigner functions 
from the formal solution of the energy stargenvalue equation (3) for the simple harmonic 
oscillator. 

The Hamiltonian reads H = \(p 2 +q 2 ), where to make it simpler we made m — h — w — 1. 
Consequently: 

A = ^I 2x2 , b T = (0,0), AJA=- A J. (102) 
Since a — |, it follows immediately from (99,100) that: 

+oo 

A,(ff - £) = 2vr £ FZ(q, P )S (E - E n ) , (103) 

n=0 



where: 



FZ(q, P )= K -^e- 2H L n (4H), E n = n+^-. (104) 

7T I 

Let us now derive an expression for the non-diagonal elements F^ n (q,p). The creation 
and annihilation operators are defined by: 



q + ip „ + q — ip 



a, 



1. (105) 



Setting E' = E + n, eq.(34) reads: 

A* (if, E', E) = W(\E' >< E\) = Jdajdk 5^\a)e^ H - E ^ a \ (106) 

The evaluation of the *-exponential yields: 

gttH+o^t = £ aa^ {q-^P+^j, (107) 

where <p(q,p) = el kH . Consequently, 

A* (if, E', E) = ±- fdk e- lkE ^ in \k), (108) 

where ^ n \k) is an analytical periodic distribution (with period An) in the variable k e 
C\{(2n+l)vr, n E Z\. 



ij {n \k) 



da 

Its Fourier expansion is: 



_aat(l-itan(fc/2)) 

f e 2iHtan(fe/2) 

cos(fc/2) 



(109) 

a=0 



+oo 



^«)(A;) = ^rf_ 2 ,_ ie -^+ 1 ) fc /2 ? (no) 
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with: 



d-a- 1 (q,p) = ±f c dz ^- 1 (, 2 + 1)- 1 (^ T ) e 



( 2a t)« 



-e- 2ff |^ J^ds z -2(/+n)-i (z 2 + !)-i e ^T = (2at)"e- 2 ^2(-l) / Lr(4iJ). 

(111) 

Substituting (110) in eq.(108), we get: 



+00 



A*(H, E', E) = J2HE~ E t ) d-a-L (112) 
1=0 

The normalized Wigner function F^ h (with m = n + I and f3 m j = J -^7) is given by: 



FZ(Q,P) = §^-21-1 = — (-WA (2a t ) m " ,i r , (4^, (m > Z > 0). (113) 

/7T 7T V 777.1 \ ' 



The analogous result for / > m is tantamount to performing the substitutions / < > m and 

< > a in the previous expression. This can be summarized in the following formula: 

,,, , , -2H , > / , 1 in m / I \ ksiqn(m — n) , , \m-n\+m-n 

FL(q,P) = e—(-l)in+m-\m-n |)/2 ^ 2 (2fl) a x 

(114) 

|m — n| — rra 

x(2at) 2 L '^-'i^i (4ff), n,m = 0,l,2,... 

2 

Notice that (-P^^p)) = F^ n (q,p), which is agreement with the quantum mechanical 

counterpart (\n >< m\) = \m >< n\. A lengthy calculation shows that F^ m obeys the 
correct *-genvalue equations: 

H*FZ l = E n FZ l , F% n *H = E m F% a . (115) 



6 The linear potential 

The hamiltonian for a particle coupled to a linear potential has a continuous spectrum. It 
is defined by: 

H(q,p) = ^ + q. (116) 

We then have: 

A=(l °), & T = (0,1) (117) 
Notice that the matrix A is singular. To solve this problem, let us define: 




(118) 
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We shall compute the *-exponential for finite A and then eventually take the limit A — > 0. 
From these expressions we get: 



b T A- x 1 b = j, A X JA X = ^J 



IX 



And so: 

AkH 



= lim [cos (a\k)]~ exp <^ — tan (ct\h) H — - — [tan (a\k) — a\k] \ = 



It then follows that: 



A* (#-£) = [dk e lk{H ~ E)+i £ = 2 At [2 (H - E)] , 
2ir J 

where Ai(z) is the Airy function [35]: 



Ai(x) = — I dk e 



-ikx-^k 3 



(119) 



(120) 



(121) 



(122) 



This is a well-known result [20] . To compute the non-diagonal elements, we need the operator 
B which implements the infinitesimal translations in the spectrum of H . A straightforward 
calculation shows that B = p. We then have to compute the *-exponential (87), with: 



A 



2 



P = ik, A\ = 
After a simple calculation eq.(15) reads: 

A* (if, E', E) = 2e- l(E '- E)p Ai 



i I b 1 = --(E'-E),l . 



2 H- 







E + E r 



(123) 



(124) 



Notice that the previous expression reduces to (121) if E' = E. Moreover using the fact that 



Ai(x) satisfies the equation, j^Ai(x) = xAi(x), we can prove that: 



dx' 



j + qj * A, (H, E', E) = E'A* (if, E' , E) , A, (H, E' , E) * ^ + q^j = EA* (H, E', E) . 

(125) 



7 Two-dimensional harmonic oscillator 

The Hamiltonian for this model is: 

^ = p| p| f f (126) 

2 2 2 2' y ! 

Our purpose is to determine the simultaneous *-genfunctions of the hamiltonian and of the 
z component of the angular momentum: 

L 3 = qifa -Pifa- (127) 
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Since the two observables commute our first step is to compute the *-exponential e l ^ H + lk ^ . 
Let us define: A = A 1 + A 2 , where: 

A l = a/ 4x4 , A 2 = cV = c ^ ^ ^ . (128) 

Here, a = c = y an d — 1> 2, 3) are Pauli's spin matrices. 
Now consider the matrices: 

S, = v^/4x4, S 2 = f-( p ~{ 2 * 2 ) . (129) 

V 2 \ -<2x2 ^2 / 

They satisfy: 

£([S 1 = A 1 , S 2 r S 2 = A 2 , SfS 2 + S%S 1 = 0. (130) 
We conclude that S A = Sx + S 2 . Moreover we have: 

SiJS^ + S 2 JSf = 0. (131) 
It then follows that B — B\ + -B 2 , where: 

Bi = -aJ, 5 2 = aR = c f ~^ 2 _^ J . (132) 



Consequently: 



= (ia) 2n h x4 , B 2n+l = i(ia) 2n+1 J, 
B 2n = c 2n h xA , B 2n+1 = c 2n+1 R, 



\ R2n+1 _ „2n+l d n — 0, 1, 2, • • • (133) 



This entails that: 



sin_Bi = —2 sin (jf) J, cosB 1 = cos 0$f) J 4x4 , 
sin_B 2 = sin (Jf^ R, cos B 2 = cos (jfj J 4x4 . 



It is easy to check that B x and B 2 commute. It then follows that: 



(134) 



sin B = sin B 1 cos B 2 + cos B\ sin £? 2 = — i sin fy ) cos (jf^j J + cos (y) sin (jf^j R, 
cos 5 = cos Bi cos 5 2 - sin B x sin _B 2 = cos (Jfj cos (Jfj J 4x4 + i sin (y) sin (y) V. 

(135) 

To compute tan B, we first need (cos B) : 

(coz B) -i cos (j) cos (j) J ^ - ! sin (j) sin (|) v fl36) 

1 ' cos 2 (y) cos 2 (y ) - sin 2 (y) sin 2 (f ) ' 1 ' 

We then get: 

/ h + k 2 \ (R- iJ\ (ki-k 2 \ /R + iJ\ . „, 

tan£ = sinE (cos5) 1 = tan — - — - tan — - — . (137) 
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It remains to determine the inverse of Sa- First of all notice that: 



A ~" = T^T2 ( %k l hx2 ~^T 2 ) • ( 138 ) 
k\-k\ \ k 2 a 2 tkj 2x2 J 

Since A = S^Sa, we conclude that S^ 1 = A^S^. From this follows that JS^J = — S^ 1 . 
Collecting all the results we get from (55): 



os(^)cos(^)~ 

(139) 

x exp [i tan (H + L 3 ) + i tan (^f^) (H - L 3 ) } . 

Let us define x = fcl + fc2 , y = kl ~ k2 and (f)(x,y) = tp(ki,k 2 ). Again this is a distribution, 
periodic in the variables x, y, with periods 2tt, and analytic for x ^ (2n + l)vr/2, y ^ 
(2m + 1)tt/2 (n, m G Z). Expanding in Fourier modes, we have: 

+oo 

<j>(x,y) = £ C n , m e inx+imy . (140) 

n,m=— oo 

The Fourier modes are given by 

C n , m = — ^ f rfx f rfy e- ira - im ^(x, y) = n n (if + L 3 ) II m (H - L 3 ) , (141) 

{2,11) J-n J-n 

where 



n„(-u) = — / cfe = — i dz z~ n (z 2 + 1) exp 

2ir J-n cos a; in Jc v ' 



(142) 



We conclude that H n (u) = for n < or n even, and that: 

U 2n+1 (u) =2(-l) n e- u L n (2u), ra = 0,1,2,-.. (143) 

We then have: 

A„(lf - £) * A*(L 3 - F) = ^ J dh J dk 2 e f^-^^-F) = 
= *fdxfdy <j)(x,y)e-^ E ~^ F = 

(144) 

= 2E n % o n 2n+1 (if + L 3 )n 2m+1 (# - L 3 )5(n — E — F)5{m — E + F) — 

= 2 Et" E ae / r n r+s+ i(H + L 3 )II r _ s+1 (i7 - L 3 )5(£ - E r )5(F - F s ), 

where r = n + m, s = n — m, E r = r + 1, F s = s and I r = {r, r — 2, r — 4, • • • , 2 — r, — r}. 
It is interesting to note that had we computed A* (if — E), A*(L 3 — F) separately, then the 
two spectra would have been independent. It is the *-product between the two expressions 
which forces the spectrum of L 3 to depend upon that of H. 
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Let us now compute the non-diagonal elements. The translation operator for L 3 is 

T = v\ - v\ + ql - q 2 2 + 2 ^i^2 + 2ipip 2 , (145) 

which satisfies: 



H,T = 0, L 3 ,T 



= 2T. 



We then get from (18,34): 



A* (if, L 3 ; E, E; F', F) = <— ^ J da J dh J dk 2 5^\a)e^ 



(H-E)+ik 2 (L 3 -F)+^aT 



(146) 



, (147) 



where F' = F + 2n and fi(k 2 ) = Jit?_ v To compute the non-commutative exponential, let 
us define: 

7 ahx2 + «7(°"i - ^3) co 2 \ 

-co 2 a/ 2 x2 + 17(01 - io- 3 ) J ' 



A 



(148) 



where a = c = -y and 7 = an(k 2 ). The matrix reads: 



& = W2x2 + Ml - fy)0"2 + &20"3, 

d = -b I 2x2 ~ Mi + & O2 + b 2 (J 3 . 



where: 



Here, bi = b 2 = ^ and 6 is any solution of: 



16c 2 ^ + 8 7 2 a& 2 + 7 4 = 0. 



-co 2 a/ 2x2 + 17(0-1 + io 3 ) 

-a/ 2x2 - 27(01 + ia 3 ) -co 2 



We then get from (56): 
B = 

where we used OjOj = <%/ 2x2 + ie^k^k- Let us write: B = B\ + B 2 + B 3 , with: 
i?i = —aJ, B 2 = cR, B 3 = i'j 



(149) 



(150) 



(151) 



(152) 



Oi + ia 3 
— 01 — io 3 



(153) 



We conclude that eqs. (133,134) still hold. Notice also that B 2 B 3 + B 3 B 2 = 0, 5 2 = 0. It 
then follows that: 

e ±i(B 2 +B 3 ) = e ±iB 2 ± iBaB -l sin ( 154 ) 

where B^ x = j^R- Consequently: 



' sin(£ 2 + B 3 ) = sin (f ) (r + ^B 3 ) , 
cos(B 2 + £? 3 ) = cos£? 2 = cos (^j J 4x4 . 
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(155) 



and 



Since B\ and B 2 + B 3 commute, we get: 

sini? = sin_Bi cos(-B 2 + B 3 ) + cosi?i sm(B 2 + B 3 ) = 

= -i sin cos (f ) J + cos (f ) sin (f ) (i? + ^5 3 ) , 

cos 5 = cos Bi cos(B 2 + B 3 ) — sin .£?! sin(S 2 + B 3 ) = 
= cos (^) cos (f ) J 4X 4 + i sin ) sin (f ) (V + ^z) , 
^ / <Ji + icr 3 \ 

V <Tl+i<T3 J' 

and V was denned in eq.(128). A lengthy calculation then yields: 

(cosier 1 = cos( ^ a) 1 cos( ^ a) x {cos (!f ) COS (f ) J 4X4 

-.sin(f) sin (f) 

and: 



where: 



(156) 



(157) 



(158) 



As before: (det cos 5) 1/2 = cos (^f^) cos (^f^). Finally: 

A-^^-^f ^ " <7(^l " <^3) T ^ 2 \ £-1^-1^ (160) 

a 2 + c 2 V ca 2 0/2x2-17(^1-^3) ) 

and JS^J = S^ 1 . Substituting these results in eqs.(55,147), we get: 

A.(H,L 3 ;E,E-F',F) = ^ J dx J dy e' l ^ E -^ F ^ n \x,y), (161) 

where 

1 d n 

ip( n) (x,y) = — {i(H + L 3 ) tanx + i(H - L 3 ) tany + aTYtana; + i)(tany - i)} 

cos x cos 1/ aa" 

(162) 

is a periodic distribution in the variables x, y with periods 2n, and analytic for x ^ (21 + 
l)7r/2, y 7^ (2m + l)7r/2 (Z,m G Z). Again we expand it in a Fourier series: 

+00 

^ n \x,y)= ]T C^e ilx+im y, (163) 

l,m=— 00 

where 

C ^ = JwL dX L dy ^ n) ^yy- llX - imV = 7^^(11 + L3)J£\H-L 3 ). (164) 



a=0 
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Here, 



r(«) 



f-« dx^^ exp [i(H + L 3 ) tan* - ilx] 



= -{2i) n+1 § c dzz-\z 2 + I)"™" 1 exp [(if + L 3 ) ($=±] 
Again, this vanishes if Z < or if / is even. We are left with: 

4^ = (22) n+1 (-l)' +1 2«e-( //+L3 )L[ l (2iJ + 2L 3 ), Z = 0, 1, 2, 
Likewise, the nonvanishing terms of Jffl are given by: 



J. 



(n) 



2m+l 



(tan y—i) r ' 
cosy 



exp — L 3 ) tany — 2 (2m + l)y] 



(-2i) n+1 e" (// - L3) 27ri(-l) m - n L™ _„(2iJ - 2L 3 ), m > 



(165) 



(166) 



(167) 



n. 



We then get: 



A*(H,L 3 ;E,E;F',F) = ^T^Tl^ f dx f dy c£L ^ e -™(E + F-2i-i)-i y (E-F-2m-i) 



- 2tt 2 ^1^0 Sm=ii I dx J dy C 2 [ +lt 2m+l ( 

E^E s eJ r - E r )5(F - F s )C^ s+1>r _ s+1 , 



(168) 

where J r = {— r, — r + 2, — r + 4, • • • , t r } and t r = r — n (for n even) or r — n — 1 (for n odd). 

" = 2T and " 



Using the fact that L 3 ,T 

\-r + 2l> 
This means that (cf.(29)): 

P- 



[r-l)\ (T 



\ (r-l')\ \ A 



T, = 16L 3 , we can prove that: 
l-V 

\-r + 2l'>, l>l'. 



(169) 



r+l+l',-r+2l 



(r-l)\ n 

x 



\(r-l')\ U 



i-v 



The non-diagonal Wigner functions F™, TS can be read off from eqs. (168,169)) for s' > s. 
To obtain F^, rs for s' < s, we complex conjugate the latter expression and perform the 
substitution s' < — ► s. The general result (for arbitrary m and uj and restoring the TVs) can 
then be written in compact form: 



("l) r 



rs',rs (Trh) 2 uj 



sign(s—s')/2 



I s — s +s — s 
4 



/ t \ 4 / Tt y 



2H 

e un. x 



(170) 



|s — s f I 



yf 2 

^ 2r+s + s'-|s-s 



/ 2(g+q;L3) \ r^-" / 2(H-^L 3 ) \ 



with s, s' E I r — {— r, — r + 2, • • • , r — 2, r}. A lengthy calculation shows that: 

H * Frs',rs = Frs>,rs * H = E r F^, rs , 

< 

k L3 * F^>, rs — Fs>F™, rs , F™, rs * L 3 = F s F^, rs , 
with E r — uh (r + |) and F s = sk. 



(171) 
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8 Conclusions 



Let us briefly summarize our results. We derived a formal solution for an arbitrary (diagonal 
and non-diagonal) *-genvalue equation. The formalism was designed to incorporate observ- 
ables with discrete and continuous spectra as well as multidimensional systems. This formal 
solution was interpreted as a non-commutative generalization of the Dirac delta distribution, 
identifying *-hypersurfaces in the non-commutative phase space. Moreover, we presented the 
complete specification of Wigner quantum mechanics in the Heisenberg picture. We further 
derived the explicit expression of our formal solution for an arbitrary quadratic phase space 
functional. If the Hessian matrix of this functional is proportional to a symplectic matrix, 
then the aforementioned explicit expression simplifies drastically and we can systematically 
analyze the continuous and discrete spectrum cases. Finally, we apply our formalism to the 
one- and two-dimensional harmonic oscillators and to the linear potential. The diagonal 
terms were previously solved in the literature ([20], [23], [9]) using other methods. The 
non-diagonal terms are to the best of our knowledge new. 
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